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Abstract 

We present an exact analytic solution to the neutrino evolution equation in mat- 
ter with periodic step-function density profile and discuss in detail the parametric 
resonance of neutrino oscillations that can occur in such a system. Solar and at- 
mospheric neutrinos traversing the earth pass through layers of alternating density 
and can therefore experience parametric resonance of their oscillations. Atmospheric 
neutrinos can undergo parametrically enhanced oscillations in the earth when their 
trajectories deviate from the vertical by about 26° — 32°. Solar neutrinos travers- 
ing the earth can experience a strong parametric resonance of their oscillations in a 
wide range of zenith angles. If the small mixing angle MSW effect is the solution of 
the solar neutrino problem, the oscillations of solar neutrinos crossing the core of the 
earth must undergo strong parametric resonance; this phenomenon should facilitate 
significantly the observation of the day-night effect in oscillations of solar neutrinos. 
If observed, the enhanced day- night effect for core crossing neutrinos would therefore 
confirm both the MSW solution of the solar neutrino problem and the parametric 
resonance of neutrino oscillations. 
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1 Introduction 



It is well known that the resonantly enhanced neutrino oscillations in matter, the Mikheyev- 
Smirnov-Wolfenstein (MSW) effect [p], §, have a simple analogue in classical mechanics: 
oscillations of two weakly coupled pendulums with slowly changing frequencies [^, |J. One 
then naturally wonders if there are any other resonance phenomena in mechanics which 
might have analogues in neutrino physics. 

One such phenomenon is parametric resonance. The parametric resonance can occur 
in dynamical systems whose parameters vary periodically with time. A textbook example 
is a pendulum with vertically oscillating point of support |5|, Under certain conditions 
topmost, normally unstable, equilibrium point becomes stable. The pendulum can oscillate 
around this point in the upside-down position. Another example, familiar to everybody, 
is a swing, which is just a pendulum with periodically changing effective length. It is the 
parametric resonance that makes it possible to swing on a swing. 

What would be an analogue of the parametric resonance for neutrino systems? Since 
matter affects neutrino oscillations, periodically varying conditions can be achieved if a beam 
of oscillating neutrinos propagates through a medium with periodically modulated density. 
For certain relations between the period of density modulation and oscillation length, the 
parametric resonance occurs and the oscillations get enhanced. The probability of neutrino 
transition from one flavor state to another may become close to unity even for small mixing 
angle. 

This phenomenon is very different from the MSW effect. Indeed, at the MSW resonance 
the neutrino mixing in matter becomes maximal (8 m = vr/4) even if the vacuum mixing angle 
8 vac is small. This leads to large-amplitude neutrino oscillations in a matter of constant 
density equal (or almost equal) to the resonance density, or to a strong flavor conversion in 
the case of matter density slowly varying along the neutrino path and passing through the 
resonance value. The situation is quite different in the case of the parametric resonance. 
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The mixing angle in matter (the oscillation amplitude) does not in general become large. 
What happens is an amplification of the transition probability because of specific phase 
relationships. Thus, in the case of the parametric resonance it is the phase of oscillations 
(rather than their amplitude) that undergoes important modification. 

The parametric resonance of neutrino oscillations has a very simple physical interpre- 
tation: the average value of the transition probability, around which the oscillations occur, 
drifts. This may lead to large probabilities of flavor transitions. We shall discuss this 
interpretation in the last section of this paper. 

The possibility of the parametric resonance of neutrino oscillations was suggested inde- 
pendently in |7| and ||. In ref. an approximate solution for sinusoidal matter density 
profile was found. In || an exact analytic solution for the periodic step-function density 
profile was obtained. However, in ref. || the results were presented only for a simplified 
case of small matter effects on the oscillations length and mixing angle, l m ~ l vac , 8 m w 9 vac . 
Parametric effects in neutrino oscillations were further studied in Q where combined ac- 
tion of parametric and MSW resonances and possible consequences for solar and supernova 
neutrinos were discussed. 

Recently, there has been an increasing interest in parametric resonance of neutrino 
oscillations. It was pointed out JCTI] that atmospheric neutrinos traversing the earth travel 
through layers of alternating density and can therefore undergo parametrically enhanced 
oscillations. The same was also shown to be true for solar neutrinos passing through the 
earth |TT]. Interestingly, this situation very closely corresponds to a periodic step-function 



density profile studied in ||. We therefore believe that it would be useful to present the 
exact analytic solution for this case in full, without assuming smallness of matter effects. We 
also study the evolution of oscillating neutrinos in the earth and derive an exact analytic 
expression for the neutrino evolution matrix in the constant-density-layers model of the 
earth structure. Although in this case neutrinos travel only over "one and a half" periods 
of density modulation, parametric resonance effects are possible and can be very large in 
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this case, too. We discuss these effects and their possible manifestations for solar and 
atmospheric neutrinos. 

This paper is organized as follows. The evolution equation for a system of oscillating 
neutrinos in matter is reviewed and the solution for the periodic step-function density profile 
is derived in Sec. 2. Various limiting cases are considered and the parametric resonance of 
neutrino oscillations is studied in Sec. 3. Parametric enhancement of oscillations of solar 
and atmospheric neutrinos in the matter of the earth is considered in Sec. 4. The results 
are discussed and a simple physical interpretation of the parametric resonance in neutrino 
oscillations is given Sec. 5. 



2 Evolution equation and its solution 



Consider a system of two relativistic neutrino species with mixing. The evolution of a such 
a system in matter is described in the weak eigenstate basis by the Schrodinger equation [0 
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Here v a ,b{t) are the probability amplitudes of finding neutrinos of the corresponding flavor 
a, b at a time t (in particular, one of these two species can be a sterile neutrino v s ). The 
parameters A and B are 
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Here Gp is the Fermi constant, E is neutrino energy, Am 2 



nr. 



— m 2 , where m 12 are the 



neutrino mass eigenvalues, and 8 vac is the mixing angle in vacuum. The effective density 
N(t) depends on the type of the neutrinos taking part in the oscillations: 

N e for v e 



N 







for v.. 



N e -N n /2 for v e 
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for v t 



Here N e and N n are the electron and neutron number densities, respectively. For transitions 
between antineutrinos one should substitute —N for N in eq. (fj). If overall matter density 
and/or chemical composition varies along the neutrino path, the effective density iV depends 
on the neutrino coordinate t. The instantaneous oscillation length l m (t) and mixing angle 
9 m (t) in matter are given by 



l m (t)=Tr/u(t), am2e m (t)=B/u(t), u(t) = yj B 2 + A(t) 2 . (4) 

The MSW resonance corresponds to A(t res ) = 0, sm29 m (t res ) = 1. 

We now consider a special case when the effective density N(t) (and therefore A(t)) is 
a periodic step function: 



N(t) 



Nt, < t < Ti 

N 2 , T x < t < T x + T 2 



N(t + T)= N(t) , T = T 1 + T 2 . (5) 

Here Ni and ^2 are constants. We shall call this the "castle wall" density profile. The 
function A(t) is expressed by a similar formula with constants A\ and A 2 . Thus, the 

Hamiltonian H(t) is also a periodic function of time with the period T: H(t + T) = H(t). 
Let us denote 

s = ^- V ' = T2 N ' (i = 1,2) - (6) 

In this notation 



A { = cos 26 vac 5 - Vi , B = sin 28 vac 5 , u { = ^(cos28 vac 5 - VA 2 + (sin 29 vac 5) 2 . 

(7) 

Any instant of time in the evolution of the neutrino system belongs to one of the two 
kinds of the time intervals: 

(1) : + nT<t<T 1 + nT 

(8) 

(2) : Ti + nT <t <T 1 +T 2 + nT , n = 0, 1, 2, ... 
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In either of these time intervals the Hamiltonian H is a constant matrix which we denote 
H\ and H 2 , respectively. Let us define the evolution matrices for the intervals of time (0, Tj.) 
and (T 1 ,T 1 + T 2 ): 

U x = exp(-i J ff 1 T 1 ) , U 2 = expH# 2 T 2 ) (9) 
The evolution matrix for the period is then 

U T = U 2 U X . (10) 

It is convenient to rewrite the Hamiltonian H given by eq. (|l|) using Pauli's a matrices: 

H{t) = Ba x - A(t)a 3 , (11) 

Let us introduce the unit vectors 

nx = —(B,0, - A x ) = (sin 20i, 0, -cos 2^), 

n 2 = —(B, 0, - A 2 ) = (sin 20 2 , 0, - cos20 2 ) , (12) 

where B\ )2 are the mixing angles in matter at densities N\ and N 2 : Q\ = 8 m (N x ), 9 2 = m (N 2 ). 
Then one can write 

Hi = cjj(crnj) . (13) 

Notice that eigenvalues of Hi are Using eqs. (|9|) and fllCf ) one arrives at the following 
expression for the matrix of evolution over the period: 

U T = Y - iaX = exp[-i(<rX)$] . (14) 

Here 

Y = cxc 2 - (n 1 n 2 )s 1 s 2 , (15) 

X = sic 2 ni + s 2 Ci n 2 - SiS 2 (ni x n 2 ) , (16) 

X 

$ = arccosF = arcsinX , X = — , (17) 
and we have used the notation 

Si = sin 4>i, d = cos 4>i , 4>i = UiTi (i = 1, 2) . (18) 



Notice that Y 2 + X 2 = 1 as a consequence of unitarity of Ut- From eq. (0) one has 

nm 2 = -^—{B 2 + A X A 2 ) = cos(20i - 26 2 ) , (19) 

(m x n 2 ) = — *— (0, B{A 2 - At), 0) = (0, sin(2^ - 20 2 ), 0) . (20) 
The vector X can be written in components as 

X= LB 1 , (Ai-A 2 ), - — sic 2 H s 2 ci . (21) 

V V iO\ 10% / Wio; 2 \a;i u; 2 // 

The evolution matrix for n periods (n=l, 2,...) can be obtained by raising Ut to the 
nth power: 

U nT = U(t = nT, 0) = exp[-i(crX)n$] . (22) 
Eqs. (|l4])-(|2~2|) give the exact solution of the evolution equation for any instant of time that 



is an integer multiple of the period T. In order to obtain the solution for nT < t < (n+ 1)T 
one has to evolve the solution at t — nT by applying the evolution matrix 

C/i(t, nT) = exp[-iffi ■ (t - nT)) (23) 

for nT < t < nT + T\ or 

C/ 2 (t, nT + T 1 )U 1 = exp[-iH 2 ■ (t - nT - T x )] expl-iH^] (24) 
for nT + T\ < t < (n + 1)T, with Hi 2 given by eq. (O). 



In the limit v4 2 = Ai eqs. reproduce the well known results for neutrino oscil- 

lations in a medium of uniform density, as they should. 



3 Parametric resonance 

Assume that the initial neutrino state at t = is u a . The probability of finding v\> at a time 
t > (transition probability) is then P{y a — > 1/5, t) = |f/ 2 i(t)| 2 where f/(t) is the evolution 
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matrix. We will concentrate now on neutrino transition probabilities for t = nT. Having 
found it, one can apply the procedure described in the end of the previous section to find 
P(v a —>■ Vb, t) for an arbitrary t. 

From eq. ( p2|) one finds 

P(y a ^u b ,t = nT) = x *lx*lxt = X^X^Xl ^ ' (25) 

This expression resembles the neutrino oscillation probability in a matter of constant density. 
However, the pre-sine factor and the argument of the sine in ( p5| ) are different from the 
amplitude and phase of neutrino oscillations in matter of either of the two densities, N\ or 
N2. In particular, the pre-sine factor need not be small even when both sin 2 26\ and sin 2 20 2 
are small. This is a consequence of the parametric enhancement of neutrino oscillations. 
The parametric resonance occurs when the pre-sine factor in (|25f) becomes equal to unity. 
The resonance condition is therefore 

Xj= ( ^ Sl c 2 + —s 2Cl ) =0. (26) 



For given values of Ni, N 2 , T\ and T 2 this equation determines the resonance energy. At 
the resonance, X 2 takes the value 



2 22 B s i 

{X ) res = (X 1 + X 2 )res = ~ — 2 (Ai — A 2 ) 

UJ 1 i0 2 



2B 2 



2„2 

(27) 

2 



! + C 2^2 



Eqs. (|25| ) -(|27 D give a general description of the parametric resonance of neutrino oscillations 
for the "castle wall" density profile of eq. (|5|) . We shall now analyze the resonance in several 
particular cases. 

3.1 Large \A t \ 

Assume first that the oscillation frequencies u)\^ are dominated by the A\^ 2 terms, i.e. 

A 2 >5 2 , (28) 
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This condition is satisfied in the following cases: 

(i) small vacuum mixing angle and the system is far from the MSW resonance, or 

(ii) (Gf/ V^)Ni : 2 ^> Am 2 /4£. The condition ( p8|) is satisfied in this case irrespective 
of the value of 9 vac . 

Eq. (^) ensures smallness of the oscillation amplitudes in matter with effective densities 
N% and N 2 . We shall consider two distinct possibilities now. 

(A) Same sign A\ and A 2 . 

This corresponds to Ni, N 2 < N MSW or Ni, N 2 > N MSW where N MS w — cos29 vac 5(\/2/GF) 
is the MSW resonance density. In this case Xf ~ sin 2 (0i+02), and the parametric resonance 
condition (E6|) becomes 



<p 1 + <j> 2 = OT = hr, k = l,2,... (29) 
Here uj is the mean oscillation frequency, 

UJ = + U 2 j; . (30) 



One can also write the resonance condition (29) as 



n = f, n^. pi) 

This is nothing but the well known parametric resonance condition in the case of small- 
amplitude oscillations |5|, ||]. Eq. (p9[) can also be written as 

T - 71 - T 

I'm — ~ — ~T i \y^) 
uj k 

i.e. the parametric resonance occurs when the mean oscillation length in matter l m equals 
the period of the density modulation divided by an integer. 

(B) Opposite sign Ai and A 2 . 

This corresponds to N\ < N MS w < ^2 or N 2 < N MS w < Ni. In this case one gets 
Xl ~ sin 2 (0! — 4> 2 ), and the resonance condition is 

0! - 2 = k'-K , fc' = 0,±l,... (33) 
8 



or 

AT Au\ k'Q 

w °-r + — ) = — • (34) 

Here we have used 

0i - 02 = AT + ^ Auj T , uj = Ul ^ ^ 2 , AT = T\ — T 2 , Au = lo\ — u 2 . (35) 



For both cases (A) and (B), at the resonance one has 



(X 2 ) res = (X 2 + X 2 2 ) res = {A l - A 2 f = sj sin 2 (2^ - 29 2 ) . (36) 

io 1 u> 2 



This follows from eqs. (|27|) and ([^). In the vicinity of the resonance sin 2 (0! ± q 
[(0i ± 02) — kit] 2 , and the transition probability can be written as 

P{v. t = «T) K -,f* „ (v^l + - t) . (37) 



Here 



D\ = -L si sin 2 (2^ - 2fl 2 ) ~ s 2 sin 2 (2^ - 2£ 2 ) , (38) 

Q = {ujiTi + uj 2 T 2 )/T , A: = 1,2, ... same sign A 1 and A 2 
u e = { , (39) 

co> AT/T + Au/2 , = 0, ±1, ... opposite sign Ai and A 2 

and we have used the fact that in the case under consideration X -C 1, so that $ ~ X. 



The conditions (|29| ) and ( |33"D determine the sum or the difference of the phases 0i and 
2 in the cases (A) and (B) respectively, but not the values of the phases themselves. These 
values can be fixed through the following consideration. The parametric resonance condition 
ensures that the pre-sine factor in eqs. (|25j ) and fl3*T|) is equal to unity. This, however, is 
not sufficient for the transition probability to be appreciable. For this one should also 
require that the argument of the sine be not too small. At the resonance, this argument is 
proportional to D^, so the fastest growth of the transition probability is realized when \D^\ 
reaches its maximum. As follows from (|3"8"|), for fixed values of mixing angles in matter 9\ 2 



this amounts to s\ — ±1, i.e. 

0! = tt/2 + fcV . (40) 



Thus, the optimal conditions for having a sizeable transition probability are realized when 
there is a parametric resonance and in addition D\ reaches its maximum. This can be 
summarized as 

1 = | + A;7r, 2 = | + fc'vr , k, k' = 0, 1, 2, ... (41) 



which follows from eqs. fl33|) and (|4"0|) . The conditions (|4]) apply to both cases (A) 

and (B). 

The parameter also determines the width of the resonance. The pre-sine factor 

in ( P?D becomes equal to 1/2 when the detuning (u e — kQ/2) = ±-Dfc. If the parametric 

resonance condition is considered as a condition on neutrino energy, determines the 

energy width of the resonance. It is easy to find the energy width at half height AE. The 

condition ( P%D is satisfied provided that either sin 2 29 vac <C 1 (barring the proximity to the 

MSW resonance), or Vi ^> 5. In the first well as in the second case for not too large 

vacuum mixing angles (cos 29 vac > 5/2Vi <C 1) one finds 

AE ^ 2D k = |2ai sin(2fli-2fl 2 )| 
E cos 28 vac S cos 26 vac 5 T 

Here Eq is the resonance energy and 5q = 5(Eq). In the case V% ^> 5q and vacuum mixing 

close to maximal the result is 

Ag^2A T Vl V 2 ^ m 

E ~ 5 5 (TiVi + T 2 V 2 ) ' ZUvac ~ 2V, K< ' ^ 

The resonance widths AE / Eq in eqs. ( f42"l) and ( |43l) are maximal when s\ = ±1, i.e. eq. 
(|40| ) is satisfied. The widths rapidly decrease with increasing resonance order fc. For the 
principal resonance, k = 1 (k = in the case of the opposite sign A\ and A 2 ), 

S -§- ~ | ffln(20! - 2# 2 ) | , 6 >Vi; ~ | sin(2^ - 2# 2 ) | H , 8 « ^ . (44) 

Consider now the transition probabilities in more detail. If the condition ( f£0| ) is satisfied, 
X res = | -Dfc | max = | sin(26 l i — 29 2 )\, $ res = 2(6*i — 6*2), and the transition probability over one 
period of density modulation T is 

P(u a -> ^, T) = sin 2 (20 a - 2fl 2 ) . (45) 
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The smallness of neutrino mixing in matter, ensured by eq. fl28|), implies that 9\ and 82 
are close to or ir/2, depending on the signs of A\ and A2. We first consider the case 
(A), when A\ and A2 are of the same sign. For A%,A 2 > (8i }2 close to 0) the probability 
(|45| ) is smaller than the maximal transition probability in matter of constant density equal 
to either Aq or N 2 , namely, sin 2 26\ or sin 2 26 l 2 . When Ai,A 2 < (8 12 close to it/ 2) the 
probability (|4"5| ) is always smaller than the largest of sin 2 28\ and sin 2 262- Thus, in this 
case the transition probability over the period cannot exceed the maximal probability in 
the case of the matter of constant density. However, the parametric resonance does lead 
to an important gain. In a medium of constant density the transition probability cannot 
exceed sin 2 29 m , no matter how long the distance that neutrinos travel. On the contrary, in 
the matter with "castle wall" density profile, if the optimal parametric resonance conditions 
( |4"T| ) are satisfied, the transition probability can become large provided neutrinos travel large 
enough distance. For t = nT the transition probability takes the value 

P{y a v b , t = nT) = sin 2 [2^! - 8 2 )\ . (46) 

This probability can become quite sizeable even for small 81 and 82 provided that neutrinos 
have traveled sufficiently large distance. The number of periods neutrinos have to pass in 
order to experience a complete (or almost complete) conversion is 

n ^whkY (47) 

However, the transition probability can be appreciable even for smaller values of n. 

Consider now the transition probabilities in the case (B), opposite sign A\ and A 2 . 
Assume for definiteness A\ > 0, A 2 < 0. Since at the MSW resonance A(t) = 0, in this 
case Ni < Nmsw < N 2 . Small mixing in matter [eq. (p8|) 1 then implies that 6\ is close to 
and 82 is close to tt/2. The transition probability over one period and n periods in the 
"optimal parametric resonance" are again given by eqs. (KSI) and ([|Bp respectively. However 
in this case, for 6 2 > it /4 + 0\/2 (which is always satisfied for small mixing in matter), one 
has sin 2 (2^2 — 20i) > sin 2 26 l i, sin 2 2# 2 . This means that even for the time interval equal 
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to one period of matter density modulation the transition probability exceeds the maximal 
probabilities of oscillations in matter of constant densities Ni and N%. The effect increases 
with increasing mixing in matter. This is a very important case; as we shall see in Sec. 4, 
this parametric enhancement is further magnified in the case of neutrinos traveling over "one 
and a half" periods of density modulation, which has important implications for neutrinos 
traversing the earth. 



3.2 Small frequency variations 

We will assume here that the frequency variations are small, oui ~ uj 2 . This is the case when 
either relative variation of density is small, \Ni — N 2 \ <§C Ni, or matter effects are small, 
Vl,2 ^ cos29 vac S. However, we will not be assuming the condition (|28[ ) in this subsection. 

Let us introduce A = (A\ + A 2 )/2 and write 

A 1 =A + a, A 2 = A -a, \a\ < \A \ . (48) 

Then 

u x ~ u ( 1 + ^r) , ^-ujJl- ^f] , u = (wi + w 2 )/2 ^ J B 2 + A 2 . (49) 
Notice that the mean oscillation frequency 

u = ^(0i + 2 ) =u (1 + O(Au/uj )). (50) 



The resonance condition fl26|) in this case coincides with eq. (P9|), but eq. ((Ej) does not 
have to satisfied. 

In the vicinity of the parametric resonance uj ~ kfl/2 one can write 

XI ~ 4 sm 2 (0! + <h) - 4^ - . (51) 

UJ Q u \ I ) 

so that the probability of the oscillations becomes 

D l _:..2 (A« 



PK ^t = nT) = gi + ( ^ m/2)2 sin' ^ ^ + <* - *"/2> 2 « ) • (52) 
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The parameter D\ here is defined through 

D k = j^2rp2 (^1 + X 2)res ■ (53) 

It coincides with the one given in eq. (|3~8D in the limit Af ^> B 2 . Consider now the case 
T\ = T 2 = T/2. We have two distinct situations. 

(i) Odd k (this includes the principal resonance, k = 1). In this case 0i — 2 — 0, and 

B 2 ( A 2 -A i y 2 
oj f UJ 2 



{X 2 f res = \ \ - x ± = sin 2 (2^ - 26 2 ) , (X 1 )L«(^ 2 )L, (54) 



2 _ &_ (A - A 2 f 

k ~Aj 7r2jfc2 ' (55) 

(ii) Even k. In this case the main terms in X 2 cancel, X 2 = 0(a 3 /Aq), and X res is 
dominated by (X 1 ) res which is of the order of (a 2 / Aq). One finds 

For the principal resonance (k = 1) the resonance width at half-height is 
AE 2 \Ax-M . 1 

~ET = ~ , , , Sm 2 °vac i = = , (57) 

E ° * + ^1 - [K>M)(£o)] 2 sin 2 20„ ac 

where Vo = (Vi + V 2 )/2. The distance that neutrinos should travel in order for the oscillation 
probability to become equal or close to one is t — nT with 

n ~ - — — . 58 

In the limit Vi 2 "C 5 the results of ref. M are recovered. 



4 Evolution of oscillating neutrinos in the earth 

The earth consists of two main structures - the mantle, with density ranging from 2.7 g/cm 3 
at the surface to 5.5 g/cm 3 at the bottom, and the core, with density ranging from 9.9 to 
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12.5 g/cm 3 (see, e.g., ||12||). The electron number fraction Y e is very close to 1/2 both in 
the mantle and in the core. The radius of the earth R = 6371 km. At the border of the 
core, which has the radius r = 3486 km, there is a sharp jump of matter density from 5.5 
to 9.9 g/cm 3 , i.e. by about a factor of two. This jump is a very important feature of the 
matter density distribution in the earth. The matter densities within the mantle and within 
the core vary little compared to the density jump on the border of the core and mantle; 
therefore to a very good approximation one can consider the earth as consisting of mantle 
and core of constant densities equal to the corresponding average densities, p m ~ 4.5 g/cm 3 
and p c ~ 11.5 g/cm 3 . Comparison of neutrino transition probabilities calculated with such 
a simplified matter density profile with those calculated with actual density profile of the 



Stacey model |L2| shows a very good agreement p0| . More recent models of earth give 
very similar profiles of the matter density distribution in the earth (for a discussion see, 
e.g., fL3|). We therefore adopt such a "constant-density-layers" model of the earth density 
distribution [] . 

Interestingly, this model of matter density profile exactly corresponds to the periodic 
step-function ("castle wall") density profile studied in ref. || and in sections 2 and 3 of 
this paper. The peculiarity of the earth's density profile is that neutrinos traversing the 
earth do not pass through many periods of density modulation; at most, they can only pass 
through, loosely speaking, "one and a half" periods (this would be exactly one and a half 
periods if the distances neutrinos travel in the mantle and core were equal). However, as we 
have seen in the previous section, even in the case of one period of density modulation the 
parametric resonance can lead to a very specific amplification of the probability of neutrino 
oscillations; especially for the case of opposite sign A\ and A2 the parametric resonance 
conditions ensure the largest increase of the transition probability over one period. As we 
shall see, this effect is further enhanced when neutrinos travel over "one and a half" periods 
of density modulation. 

1 The same approximation was previously used in [ [i4[ [l5| . A different analytic approach to propagation 
of oscillating neutrinos in the earth was developed in p6| . 
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4.1 Evolution over "one and a half" periods 

Consider now the evolution matrix for neutrinos traversing the earth. We will be interested 
in the situation when neutrinos traverse the mantle, core and then again mantle; this hap- 
pens for the neutrino trajectories with the zenith angle 180° ±33.17° (nadir angle < 33.17°). 
The evolution matrix in this case is 

U 3 = U X U 2 U X = UxU T = Z- taW , (59) 

where 

Z = 2 Cl Y - c 2 , (60) 
W = 2 Sl y ni + s 2 n 2 , (61) 
and Y has been defined in (|15"|). The vector W can be written in components as 

W = ( 2 Sl Y- + s 2 - , , - ( 2 Sl Y^ + sA)). (62) 

The transition probability in this case is simply 

P(y a -> Vb) = Wl . (63) 

Eqs. (|59|)-(p3|) give the neutrino transition probability in the earth in the constant- 
density-layers approximation in general case, i.e. for arbitrary values of oscillation pa- 
rameters and neutrino energies 0. We shall now study the conditions for the parametric 
enhancement of the oscillations. For fixed values of N\, N 2 , Am 2 /4E and 8 vac these condi- 
tions give the optimal values of the phases <pi and <p 2 , i.e. of the distances T\ and T 2 . It is 

straightforward to find extrema of Wf. There are extrema of three types. 

2 The transition probability ( |63[ ) applies when the initial neutrino state at t = consists only of v a . 
In some situations, such as for solar neutrinos, the neutrino state arriving at the surface of the earth is a 
coherent superposition of v a and Uf,. This case can be easily treated using the evolution matrix (|59|). For 
example, for the probability P2e = P(^2 ~^ v e ) relevant for oscillations of solar neutrinos in the earth JlTj 
one finds P2e = sin 2 9 vac + cos 29 vac + W\ W3 sin 29 vac . Different (but equivalent) expressions for Pi e 
were derived in |l4| and |ll[] . For small 9 vac the probability Pi e practically coincides with (^3|) . 
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(1) Extrema at c\ = c 2 = 0, i.e. <f>\ = tt/2 + far, <p 2 = tt/2 + k'n. These are exactly 
the parametric resonance conditions found in Sec. 3. Thus, in the case of "one and 
a half" periods, as well as in the case of n periods, the oscillation probability reaches an 
extremum when the parametric resonance conditions are satisfied. The extremum is a 
maximum provided that either 



fi = nm 2 - L x /lo 2 - + 8) /4 < , (64) 



or 



f 2 = n in2 - u x /u 2 + Ju\/ul + 8 /4 > , (65) 



otherwise it is a saddle point. The functions f\ and f 2 are plotted in Fig. 1 along with A\ 
and A 2 for sin 2 29 vac = 0.01. For small 6 vac the conditions (Bl) and (B31) are satisfied for 



5 < V 2 (except in a small region of 5 in the vicinity of the point where A\ = 0). We shall 
discuss these conditions in more detail below. 

At the resonance the transition probability (|63D takes the value 

P{v a -> u h ) = sin 2 2(9 2 - 200 ( Ci = c 2 = 0) . (66) 

The lowest order resonance corresponds to 0i = 2 = 7r/2. The parametric resonance of this 
type is illustrated by Fig. 2. The resonance width is rather large; the transition probability 
decreases by a factor of two for 

A0 2 , 2 = |0 li2 - tt/2| ~ - . (67) 

Thus, the resonance enhancement of neutrino oscillations can occur even for quite sizeable 
detuning of the phases </>i ;2 . 

(2) Extremum at Si = c 2 = 0. It is a maximum provided A 2 > and a saddle point 
otherwise. The transition probability is 

P(y a -> Ub) = sin 2 20 2 ( Sl = c 2 = 0) . (68) 
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This case corresponds to zero transition probability in the mantle and maximal possible 
transition probability in the core. The extremum of this type is illustrated by Fig. 3. It is 
not a parametric resonance. 

(3) Maximum at cf = sf = 1/2, S2 — 0, A\ — 0. The latter condition means that the 
MSW resonance occurs at the density N\. The transition probability is equal to unity. This 
case is illustrated by Fig. 4; it requires very special conditions and we shall not discuss it 
anymore. 

The behavior of the neutrino transition probability for v e <-* r oscillations as a function 
of the phases <pi and 4>2 for various values of 5 is shown in Figs. 2-8. For large 5 for which 
A2 > the maximum of the type (2) occurs (Fig. 3). In this domain of S the conditions 
(|64l ) and ( |B~5"D are not satisfied (see Fig. 1), so the extremum at c\ = C2 = is a saddle point 
rather than maximum. Eq. ( j64|) is fulfilled in the range 9.7 x 1CT 14 eV < 5 < 2.12 x 1CT 13 eV, 
which results in a clear parametric peak at <pi = <p2 = k/2 (Fig. 2). For 5 in the range 
7.8 x 10~ 14 eV < 5 < 9.7 x 10~ 14 eV the function f\ is positive while f2 is negative, so the 
extremum at <pi = <p2 — 7r /4 is again a saddle point; however, the value of the transition 
probability at this point is close to the maximal one (Figs. 5 and 6) except in the very close 
vicinity of the value 5 = 5\ = 8.6 x 10~ 14 eV which corresponds to the MSW resonance 
at N = N\. The transition probability for 5 = 5i is shown in Fig. 4. The maxima of 
the transition probability are located at c\ = s\ = 1/2, s 2 = 0, as discussed in the case 
(3) above. For 5 < 7.8 x 10~ 14 eV the function f 2 becomes positive and the extremum at 
<pi = <p2 = ?r/4 again becomes a maximum (Fig. 7). Notice, however, that with decreasing 5 
the height of the resonance peak rapidly decreases. The shape of the transition probability 
changes especially quickly in the region 7.5 x 1CT 14 eV < (5 < 1 x 1CT 13 eV, i.e. close to 
5\. For comparison, we also plotted in Fig. 8 the transition probability for N2 = Ni\ since 
matter density is constant, no parametric resonance occurs in this case. The probability of 
u e,n,r v s oscillations exhibits a similar pattern but for 5 scaled down by a factor of two. 

As follows from (|B"Bj), whether or not the parametric resonance can lead to a significant 
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increase of the transition probability over the time equal to "one and a half" periods of 
density modulation depends on the values of the mixing angles in matter 9\ and 62 ■ The 



expression ( |66| ) for the transition probability at the resonance has been found in [T(J and 
used there to interpret a specific enhancement of the probability of atmospheric neutrino 



oscillations in the v „ <-> v s channel in terms of the parametric resonance. However, in [10 



only the case 5 <C Vi, V2 was considered. It is easy to see that, at least for not too large 
mixing in matter, effect is in fact the largest for S in the interval V\ < cos29 vac 8 < V2 
(opposite sign A\ and A 2 ), or at least close to this interval. In this range of 5 one has 
#1 = ^mantle < 7r /4, 9 2 = 9 core > 7r /4, and the expression 9 2 — 19\ is closer to 7r/4 (which 
gives maximal transition probability) than either of 9\ and #2 provided that 9 2 > 7r/4 + 9\ , 
9 2 > 3#i (small mixing in matter). Typically, the enhancement effect is even stronger 
in this situation than in the case of oscillations over one period of density modulation. 
Indeed, in that case the relevant combination of the mixing angles was #2 — #1 [see (45) ]; for 
#2 > 7r/4 + 3^i/2 the combination of the mixing angles #2 — 2#i is closer to 7r/4 than 9 2 — 9\. 



4.2 Parametric enhancement of neutrino oscillations in the earth 

In the previous subsection we discussed the behavior of the transition probability as a 
function of the phases <pi and (f>2 treating these phases as free parameters. For neutrino 
oscillations in the earth, however, the phases (f>i and <p2 are not arbitrary: they depend on 
neutrino energy and on the lengths of neutrino paths in the mantle and core. These lengths 
vary with the zenith angle of the neutrino source, but their changes are correlated and they 
cannot take arbitrary values. It is therefore not obvious whether there are any values of the 
zenith angles and neutrino parameters for which the parametric resonance conditions can 
be satisfied. As we shall see, the answer to this question is positive. 

Consider neutrinos coming to the detector from the lower hemisphere and traversing 
the mantle, core and then again mantle. As we have mentioned, this corresponds to zenith 
angles of the neutrino source B = 180° ± 33.17°. The distance T\ that neutrinos travel in 
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the mantle (each layer) and the distance T2 that they travel in the core are related to the 
zenith angle 0, earth's radius R and core radius r by 



T 1 = R\-caaQ-y/r 2 /IP-fnn 2 Q ) , T 2 = 2R ^r 2 /R 2 - sin 2 6 . (69) 

Assuming the average matter densities in the mantle and core 4.5 and 11.5 g/cm 2 respec- 
tively, one finds for the u e <-» u^ T oscillations involving only active neutrinos the following 
values of V\ and V-i- 

Vi = 8.58 x 1(T 14 eV , V 2 = 2.19 x 1(T 13 eV . (70) 

For transitions involving sterile neutrinos v e v s and v^ T <-> u s , these parameters are a 
factor of two smaller. 

Consider the lowest-order parametric resonance, 

7T 7T 

UiT x = - , u 2 T 2 = - . (71) 

Given a value of S, these conditions yield magnitudes of T\ and T2 and so fix the zenith 
angle for which they are satisfied (provided, of course, that such an angle exists). The 
problem therefore reduces to finding values of 6 and 9 vac for which eqs. ( |7lD have solutions. 
A nontrivial point here is that both conditions in eq. ( ffTD have to be satisfied at the same 
value of G. Notice that these conditions constrain the allowed values of the vacuum mixing 
angle. Assuming that the resonance condition can be satisfied for all neutrino trajectories 
that cross the core (including the vertical ones) one can derive from the second equality in 
flTT| ) the following upper limit |: 

sin 2 2g mc < f ^ 0-04. (72) 

v 2 J max " 2 

If one excludes the zenith angles close to 180°, the constraint becomes less stringent. For 
example, for sin 2 > 0.12 one obtains sin 2 28 vac < 0.07. 



3 An analogous upper bound following from the first equality is less restrictive. 
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The analysis shows that the parametric resonance conditions for neutrino oscillations 
in the earth are indeed satisfied in a number of cases; those include v e «-> v^ T as well 
as v e ,ii,T *-> v s oscillations. A systematic study of all the cases fully covering the space 
of parameters 5, sin 2 29 vac and is beyond the scope of the present paper; here we will 
discuss only a few cases of particular interest. In what follows, we shall discuss the neutrino 
trajectories in the earth in terms of the nadir angle G n = 180° — O. 

(i) v e <-> v^ T oscillations, sin 2 26 vac <C 1 (this case is especially relevant for oscillations 
of solar neutrinos in the earth). There is a broad range of nadir angles for which the 
parametric resonance conditions (ffl[) are approximately satisfied: < @ ra < 25°. The 
maximum of transition probability corresponds to 5 ~ (1.7 — 1.9) x 10~ 13 eV, depending 
on O n . Fig. 9 shows the transition probability fl63|) as a function of 5 for sin 2 29 vac = 0.01 
and G n = 11.5°. There is a parametric resonance peak at 5 ~ 1.9 x 10~ 13 eV. A "shoulder" 
on the right slope of the peak is due to the MS W resonance in the core. The transition 
probability at the peak exceeds the biggest of sin 2 26\ and sin 2 262 (which are the maximal 
oscillation probabilities in matter of constant densities equal to Aq and N 2 ) by about a 
factor of 3. The values of 0i — tt/2 and 02 — 7r/2 for this case are shown in Fig. 10. One can 
see that these functions vanish at very close values of S, i.e there is a range of S (around the 
resonance value 80 — 1.9 x 10~ 13 eV) where both conditions in eq. (|n]) are satisfied to a 
very good accuracy. This pattern does not change in the whole interval < B n $ 25°: with 
varying B n the lowest intersection point of the curves 0i —tt/2 and 02 — tt/2 moves slightly 
along the 5 axis but remains close to it. For comparison, we have also plotted in Fig. 9 the 
transition probability for V2 = V\ which exhibits no parametric resonance. 

In fact, a strong parametric enhancement of the transition probability persists up to the 
nadir angle Q n = 32.9° (after which the resonance peak merges with the MSW resonance at 
N = Aq). The corresponding resonance values of S are 5q ~ (1.1 — 1.5) x 10~ 13 eV. In this 
region of the nadir angles, the value of the transition probability at the peak is typically 
about a factor of 1.5 - 2 smaller than it is for < n < 25°. This is related to the fact 
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that the conditions (|7T|) are only approximately satisfied: the detuning of 0i and 02 varies 
between —0.3 and —0.9 in this regime. However, since in this case the resonance energy 
is farther away from the energies corresponding to the MSW resonance in the mantle and 
core, the relative enhancement of the transition probability is even stronger: for example, 
for n = 30° the transition probability at the peak exceeds the biggest of sin 2 2#i, sin 2 29 2 
by a factor of 5. 

(ii) v e <->■ z/ MjT oscillations, sin 2 28 vac = 0(1). This case may be relevant for the large 
mixing angle MSW solution of the solar neutrino problem and is also of interest for atmo- 
spheric neutrino oscillations. For sin 2 29 vac ~ 0.6 — 0.8 there is a parametric resonance in the 
narrow range of the nadir angles G n ~ 31.1° — 32.7°. The corresponding resonance values 
of S are Sq ~ (6 — 7) x 10~ 14 eV (see Figs. 11 and 12). In Fig. 11 a narrower high peak at 
5 ~ 3.1 x 10~ 13 eV is due to the MSW effect in the core. 

(in) v e ,fi,T v s oscillations, sin 2 29 vac <C 1. At small nadir angles, Q n ^ 18°, the 
conditions (|7T|) are satisfied for 5 ~ (1.5 — 1.7) x 10~ 13 eV. However, since 5 > V 2 / cos29 vac , 
eqs. (0) and (|65|) are not fulfilled; therefore this solution corresponds to a saddle point 
rather than to a maximum. There is no parametric enhancement for this range of 5. 

For a wide range of the nadir angles, < Q n < 32.6°, there is a maximum of transition 
probability with 5 varying from 1 x 10~ 13 eV (at 6 n = 0) to 6.8 x 10~ 14 eV (at Q n = 
32.6°). The corresponding values of 0i and 02, though, are X ~ 2 ~ 7r/4 rather than 
7r/2. Nevertheless, this peak can be interpreted as a parametric resonance, at least for 
Sq not very close to the ends of the allowed interval. The point is that the width of the 
parametric resonance in the (0i, 02) plane is rather large, and detuning ~ 7r/4 just decreases 
the transition probability by about a factor of two (see eq. (|J7]) and Fig. 2). Also, the 
approximate equality of 0i and 02 indicates the parametric origin of the peak. On the other 
hand, close to the end points of the allowed intervals of 5q, i.e. 5q ~ 6.8 x 10~ 14 eV and 
5q ~ 1 x 10~ 13 eV, the peak values of the transition probabilities nearly coincide with either 
sin 2 29\ or sin 2 2#2- In the vicinity of these points an interplay of the parametric and MSW 
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resonances takes place. 

The parametric resonance for G n ~ 30.9° (5q = 7.8 x 10~ 14 eV) is illustrated by Figs. 
13 and 14. The transition probability at the peak exceeds sin 2 29\ ~ sin 2 2# 2 by about a 
factor of 2.7. 

(iv) v e,n,r v s oscillations, sin 29 vac = 1- There is a parametric resonance peak for 
n = 25.8° — 32.4° and the resonance values So = (5 — 7) x 10~ 14 eV. This is the parametric 
resonance found in HIT) . 



5 Discussion 

Neutrino oscillations can be parametrically enhanced in a medium with periodically varying 
density. The periodic step-function ("castle wall") density profile is a very simple example; 
it allows for an exact solution of the neutrino evolution equation and therefore is very illumi- 
nating. In addition, this example is of practical importance since to a good approximation 
the density profile of the earth can be considered as a step-like function with nearly constant 
densities of the mantle and core. 

The parametric resonance realizes very special conditions for neutrino oscillations, lead- 
ing to a possibility of a striking increase of the transition probabilities. Typically, neutrinos 
have to travel over many periods of density modulation in order for the parametric en- 
hancement to manifest itself. However, under certain conditions even for neutrino evolution 
during one period the parametric effects can be quite sizeable. Of course, in this case the 
neutrinos are not exposed to a periodic potential, and so one may wonder how the paramet- 
ric enhancement of neutrino oscillations can occur. This can be explained as follows. The 
parametric resonance implies that the changes of the oscillation phase and matter density 
profile with the coordinate along the neutrino path are correlated in a very special way. 
This "synchronization" allows the transition probability to overbuild after every half-period 
of density modulation (see below); if the parameters A 1 and A 2 defined in eq. (|7|) are of 
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opposite sign, and in addition the oscillation amplitudes at the densities N± and iV 2 are 
not too small, a considerable increase of the transition probability is possible even for one 
period. The parametric enhancement for one and a half periods of density modulation is 
then even larger. 

It is instructive to examine under what conditions a complete neutrino flavor conversion 
over one period of density modulation is possible. Though rather contrived, this case clarifies 
the essence of the parametric resonance of neutrino oscillations. From eq. ( |4"5| ) we find that 
a complete conversion over period T would require 

2(0! - 2 ) = ±tt/2 (73) 

in addition to the conditions (f7T|). It is easy to see that, when eq. ( [73] ) is satisfied, the 
transition probability (|45|) is equal to the sum of the amplitudes of the oscillations in matter 
of constant densities N\ and iV 2 , i.e. 

sin 2 (20i - 20 2 ) = sin 2 26 1 + sin 2 29 2 . (74) 

Since at the resonance 0! = 2 = 7r/2, we have the following pattern of neutrino oscillations 
in this case. During the first part of the period, < t < Ti, usual oscillations in matter 
of constant density N\ occur. At the time t = Ti maximal possible in this case transition 
probability is achieved, which is sin 2 20i. If the density remained constant, the transition 
probability P would start decreasing after that and would return to zero at t = 1T\. How- 
ever, the density changes to N 2 and, as a consequence of the parametric resonance, the 
transition probability continues growing instead of starting decreasing. From eqs. fl7lj) and 
(|73| ) it follows that in the time interval (Ti, Ti + T 2 ) the transition probability P undergoes 
one more half-period increase (with the amplitude sin 2 20 2 ), but starting from the initial 
value sin 2 29i and not from zero. Thus, in this idealized case the parametric resonance 
places one half-wave piece of the transition probability on the top of the other: the prob- 
ability never decreases until it reaches the maximal value equal to unity. If the condition 
( |7T| ) is only approximately satisfied, there is some decrease of the transition probability after 
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the first "half-period", but P does not reach zero, and the decrease is followed by another 
increase. What happens is essentially that the average value around which the transition 
probability oscillates starts drifting towards larger values. In this way P can become quite 
large even if the amplitude of the oscillations around the average value is small. This is 
illustrated by Fig. 15, where the dependence of the transition probability on the coordinate 
along the neutrino trajectory is shown. 

Neutrino oscillations in the earth can undergo a strong parametric enhancement. Neu- 
trinos coming to the detector from the lower hemisphere from a source with the zenith angle 
in the interval = 180° ± 33.17° traverse the earth's mantle, core and then again mantle 
and so pass "one and a half" periods of density modulation. The possibility of parametric 
enhancement of atmospheric neutrino oscillations in the earth was pointed out in [|K| where 
the case 5 = Am 2 /AE -C Vi, Vi and the oscillations in the <-> v s channel were discussed. 
It was shown that the parametric resonance can modify the zenith angle distribution of the 
atmospheric neutrino events. As follows from our consideration, the parametric resonance 
can occur for oscillations of atmospheric neutrinos in both v e and <-> v s channels 

for large enough nadir angles, Q n ~ 26° — 32°, 

The energy width of the parametric resonance AE/E is typically ~ 2 — 3. Since 
the atmospheric neutrino anomaly has been observed in a rather wide range of neutrino 
energies (sub-GeV and multi-GeV), it seems unlikely that the parametric resonance can alter 
considerably the gross features of the atmospheric neutrino oscillations, such as the allowed 
values of neutrino parameters Am 2 and sin 2 29 vac . In any case, numerical analyses of the 
atmospheric neutrino data should have automatically taken the parametric enhancement 
into account. Still, we believe that it may be worth re-analyzing the data paying more 
attention to the parametric resonance and looking for its possible manifestations. It should 
be also stressed that the parametric resonance can occur in oscillations of atmospheric 
neutrinos even in those channels that are not responsible for the atmospheric neutrino 
anomaly; such effects are potentially observable and of considerable interest. 
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The parametric resonance may be very important for oscillations of the solar neutrinos 
in the earth. For small sin 2 26 vac , in a wide range of zenith angles almost completely covering 
the earth's core, the u e <-> u^ T oscillations exhibit a strong parametric resonance with the 
peak at 5 ~ (1.7 — 1.9) x 1CT 13 eV. For typical energy of 8 B solar neutrinos E ~ 8 MeV these 
values of S give Am 2 ~ (5.5 — 6.1) x 10~ 6 eV 2 . Amazingly, this almost exactly corresponds 
to the center of the allowed interval of Am 2 for the so called "small mixing angle" MSW 
solution of the solar neutrino problem (for recent discussions see, e.g., |L8|, [13], ^9|). Strong 



parametric enhancement of the probability of the oscillations of solar neutrinos in the earth, 
by up to a factor of 7, occurs in the whole allowed region of Am 2 and sin 2 29 vac . The small 
mixing angle MSW solution gives the best fit of the available solar neutrino data and so is 
the most likely solution of the solar neutrino problem. The parametric resonance can also 
occur in the v e <-+ v s oscillations in the earth, but the effect is smaller in this case. The 
resonance values 5 ~ (6.8 — 10) x 1CT 14 eV correspond to Am 2 ~ (2.2 — 3.2) x 10~ 6 eV 2 (for 
E = 8 MeV) . These values are close to the lower bound of the allowed interval of Am 2 for 
the small mixing angle solution in this channel of oscillations. For large mixing angle the 
parametric resonance is also possible, but only in a rather limited interval of zenith angles 
of neutrino trajectories. In addition, the resonance values of Am 2 in this case are below the 
allowed interval. 

In both cases of v e <-> i/ MjT and v e «-> v s oscillations the parametric effects are strongest 
at neutrino energies which are between the energies corresponding to the MSW resonances 
in the core and in the mantle of the earth. This is illustrated by Fig. 16; one can see that 
the parametric peak at 6 = 1.5 x 10~ 13 eV is between the maxima of sin 2 26\ and sin 2 29 2 
which are due to the MSW resonances in the mantle and in the core respectively. In other 
words, the MSW resonance density for neutrinos with S ~ 1.5 X 10~ 13 eV is about 7.8 g/cm 3 ; 
this density is between the core density and the mantle density. The value of the transition 
probability at the parametric peak exceeds the amplitudes of neutrino oscillations in the 
mantle sin 2 2Q\ and in the core sin 2 2Q 2 by more than a factor of six. 
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The enhancement of the transition probability in matter in the case when the minimal 
and maximal densities satisfy Ni < Nmsw < N 2 would be of no surprise if the periodic 
density profile was, e.g., of the sinusoidal type. Then somewhere the matter density would 
be equal to the MSW resonance density, and one would expect an increase of the transition 
probability ||. However, in the case of the "castle wall" density distribution the MSW 
resonance condition is not satisfied anywhere because the density profile is discontinuous: 
the only allowed values of density are N\ and N 2 . It is quite remarkable that in this case, 
too, there is a sizeable enhancement of the transition probability compared to the case of 
matter of constant density. This effect cannot be explained by the usual MSW-type increase 
of the oscillation amplitude; it is related to a specific phase relationships which are due to 
the parametric resonance. 

One of the most promising ways of testing the MSW solution of the solar neutrino 
problem is to look for the day-night effect due to regeneration of solar u^s in the matter of 
the earth during the night [O, [21], |2T| . For the most probable small mixing angle solution, 



however, this effect is expected to be rather small - at the level of (1 — 3)%. The parametric 
resonance must increase significantly the day-night effect for neutrinos crossing the earth's 
core. Unfortunately, due to their geographical location, the existing solar neutrino detectors 
have a relatively low time during which solar neutrinos pass through the core of the earth 
to reach the detector every calendar year. The Super-Kamiokande detector has a largest 
fractional core coverage time equal to 7%. In it was suggested to build a new detector 
close to the equator in order to increase the sensitivity to the earth regeneration effect; this 
would also maximize the parametric resonance effects in oscillations of solar neutrinos in 
the earth. 

The parametric enhancement of oscillations of solar neutrinos in the earth for neutrinos 
that travel significant distances in the earth's core was discovered numerically (but not 
recognized as the parametric resonance) in [^l], |l^, [22], 18]. The authors of |18|] found that 



the day-night effect in the oscillations of solar neutrinos in the earth can be enhanced by 
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up to a factor of six when solar neutrinos cross the earth's core. These results were quite 
surprising as the largest effect resulted for neutrinos of energies which correspond to the 
MSW resonance at the densities of about (6 — 8) g/cm 3 , the values lying between the core 
and mantle densities. For neutrinos of those energies the MSW resonance is inoperative []. 
The possibility of explaining these findings in terms of the parametric resonance was pointed 



seems 



out in ]nj] , where also a different interpretation was suggested which the author of |TT[ 
to prefer. Our results support the parametric-resonance interpretation. 

It is remarkable that the parametric resonance conditions can be fulfilled for neutrino 
oscillations in the earth. An important observation here is 10, |23|] that the potentials V% 
and Vi corresponding to the matter densities in the mantle and core, inverse radius of the 
earth R~ l , and typical values of S = Am 2 /AE of interest for solar (and atmospheric) neutri- 
nos, are all of the same order of magnitude - (3 x 10 -14 - 3 x 10 -13 ) eV. It is this surprising 
coincidence that makes appreciable earth effects on oscillations of solar and atmospheric 
neutrinos possible. Even more remarkable fact is that the very specific conditions on the 
phases of neutrino oscillations in the mantle and core of the earth, <pi ~ </> 2 ~ 7r/2, which 
lead to the parametric resonance of neutrino oscillations, are also satisfied in many cases of 
interest. Those include atmospheric neutrinos, and, most importantly, solar neutrinos. 

For small vacuum mixing angle, the parametric resonance occurs for the values of 
5 = Am 2 /4:E which exactly correspond to the allowed values for the small mixing angle 
MSW effect in the sun. If the small mixing angle MSW effect is the solution of the solar 
neutrino problem, the oscillations of solar neutrinos in the earth therefore must undergo 
strong parametric resonance. The parametric resonance provides a natural interpretation 
for the increase of the night-time regeneration effect in the case of core-crossing neutrinos 
found numerically in pi], [15|, |22], By selecting only those neutrinos, one can improve 
the prospects of observing the day-night effect (though at the expense of reduced statistics). 

4 Strictly speaking, the density profile of the earth is not discontinues; the densities (6 — 8) g/cm? can be 
achieved on the border of the mantle and core where the density jumps from 5.5 to 9.9 g/cm 3 . However, 
the jump occurs in a very thin layer and so the MSW resonance does not develop. 
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If observed, the enhanced day-night effect for core crossing neutrinos would therefore con- 
firm both the M SW solution of the solar neutrino problem and the parametric resonance of 
neutrino oscillations. 

The author is grateful to Q.Y. Liu and A.Yu. Smirnov for useful discussions and to S.T. 



Petcov for informing him about the results of the paper |TT] prior to its publication 
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Figure captions 



Fig. 1. Functions fi(6) (upper solid curve) and / 2 (5) (lower solid curve) [eqs. (|5^ ) and 
(|55|)1. Also shown are Ai(5) x 10 13 (dashed line) and A2 (5) X 10 13 (short-dashed line), where 
Ai >2 are defined in eq. ©. sin 2 2# mc = 0.01, n = 11.5°. 

Fig. 2. Transition probability ( |63"D for i/ e <-> z/ MjT oscillations over "one and a half periods" 
of density modulation (t = 2T\ + T 2 ) vs oscillation phases <pi and 2 acquired during the 
time intervals T\ and T 2 . sin 2 20 MC = 0.01, 5 = 1.2 x 10~ 13 eV. 

Fig. 3. Same as Fig. 2, but for 5 = 2.5 x 10~ 13 eV. 

Fig. 4. Same as Fig. 2, but for 5 = 8.6 x 10~ 14 eV. 

Fig. 5. Transition probability ( |63|) for i/ e «-> z/ M r oscillations in the earth for core-crossing 
neutrinos vs 5 (solid curve). Also shown is the probability in the case V2 = V\ = 8.58 x 10~ 14 
eV (dashed curve). sin 2 2# mc = 0.01, Q n = 11.5°. 

Fig. 6. Phases 4>i ~ 71 /2 (dashed curve) and </> 2 — n/2 (solid curve) vs 6 for v e <-> v^ T 
oscillations in the earth, sin 2 29 vac = 0.01, 0„ = 11.5°. 

Fig. 7. Transition probability fl63| ) for v e <-» v^ iT oscillations in the earth for core-crossing 
neutrinos vs 5. sin 2 29 vac = 0.6, Q n = 32.2°. 

Fig. 8. Phases <pi — n/2 (solid curve) and 2 — it/ 2 (dashed curve) vs 5 for v e <-> v^ T 
oscillations in the earth, sin 2 26 vac = 0.6, O n = 32.2°. 

Fig. 9. Transition probability ( |6"3"D for ^e iAtiT <-> v s oscillations in the earth for core-crossing 
neutrinos vs 5. sin 2 20 vac = 0.01, 9„ = 30.9°. 

Fig. 10. Phases <fti — it/A (solid curve) and </> 2 — 7r/4 (dashed curve) vs 5 for ^ e)M , r ^ v s 
oscillations in the earth, sin 2 26 vac = 0.01, G n = 30.9°. 

Fig. 11. Transition probability ( |63| ) for z/ e «-> v^ T oscillations in the earth as a function 
of the distance t (measured in units of the earth's radius R) along the neutrino trajectory. 
5 = 1.7 x KT 13 eV, sin 2 2# mc = 0.01, Q n = 11.5°. 
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Fig. 12. Transition probability (|63"D for u e <-> v^ T oscillations in the earth, (solid curve). 
Also shown are sin 2 26 1 (dashed curve) and sin 2 2# 2 (short-dashed curve), sin 2 26 vac = 0.01, 
9 n = 28.6°. 
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Fig. 6. 
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Fig. 10. 
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